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Recently, we proposed a method for calculating the "radiation reaction" self-force exerted on a 
charged particle moving in a strong field orbit in a black hole spacetime. In this approach, one first 
calculates the contribution to the "tail" part of the self force due to each multipole mode of the 
particle's self field. A certain analytic procedure is then applied to regularize the (otherwise diver- 
gent) sum over modes. This involves the derivation of certain regularization parameters using local 
analysis of the (retarded) Green's function. In the present paper we present a detailed formulation 
of this mode-sum regularization scheme for a scalar charge on a class of static spherically-symmetric 
backgrounds (including, e.g., the Schwarzschild, Reissner-Nordstrom, and Schwarzschild-de Sitter 
spacetimes). We fully implement the regularization scheme for an arbitrary radial trajectory (not 
necessarily geodesic) by explicitly calculating all necessary regularization parameters in this case. 

(N 

04.25-g, 04.30.Db, 04.70.Bw 



00 



in 
o 
o 



I. INTRODUCTION 



The motion of a test pointlike mass (a "particle") in orbit outside a black hole is commonly studied to model, and 
gain understanding of, realistic astrophysical scenarios involving highly relativistic two-body systems — particularly, 
the capture of a small compact object by a supermassive black hole jlj]. To describe the orbital evolution of such a 
particle on a strongly curved background, one must take into account non-geodesic effects caused by the interaction 
of the particle with its own gravitational field. This problem of deducing the self force (or "radiation reaction" force) 
exerted on the particle is often treated via perturbation theory: one assumes that the particle is endowed with a 
charge q much smaller than the mass of the black hole (this charge may represent the particle's mass, electric charge 
— > • or — as in the toy model studied in the current paper — its scalar charge) and looks for the 0(q 2 ) correction to the 
equation of motion. The basic challenging task involved in this calculation, already in flat space, is, of course, correctly 
I , handling the divergence of the self field at the very location of the particle; namely, the introduction and justification 
^Jy of an appropriate regularization method. When considering the case of curved spacetime, additional difficulty arises 
. . , due to the nonlocal nature of the self-force effect: waves emitted by the particle at some moment may backscatter off 
k*" ■ spacetime curvature and interact back with the particle at later stages of its motion. The occurrence of this so-called 
"tail" contribution to the self force results in that the calculation of the self force at a given moment requires, in 
principle, knowledge of the entire causal history of the particle. A number of methods (briefly surveyed below) have 
been proposed over the years for calculating the self force in curved backgrounds. The interest in this problem has 
greatly risen lately by virtue of recent years' developments towards experimental gravitational waves detection, and 
the consequent need for accurate predictions for the orbital evolution of strongly gravitating two-body systems. Yet, 
actual calculations of the self force have been restricted thus far only to a few very simple cases (see below). 

The standard technique for calculating the radiative evolution of orbits around black holes is the one based on 
energy-momentum balance considerations S. In this approach one computes the flux to infinity (and across the 
horizon) of quantities associated with the constants of motion in the lack of self-force effects (specifically, in the 
Schwarzschild and Kerr backgrounds, the particle's energy E and azimuthal angular momentum L z ), thus deducing 
the temporal rate of change of these "constants". Such balance calculations, though developed to a great extent, 
present two basic drawbacks: First, in the important Kerr case they are inapplicable for calculating the rate of change 
of the third constant of motion necessary for a full specification of the orbital evolution, i.e., the Carter constant Q, 
as this quantity is not additive. Second, these calculations do not account for the non-dissipative, yet important, part 
of the self force §,§. 

For the above reasons, a method based on direct calculation of the momentary force along the worldline seems more 
adequate. In the context of electromagnetism in flat space, such a method is familiar from the classic work by Dirac 
H , concerning the electromagnetic self force on a (classical) electron. Dirac avoided the singularity of the self field 
at the particle's location by introducing the "radiative potential" , constructed by taking the difference between the 
retarded and advanced electromagnetic potentials (what results in the cancellation of the problematic singular part). 
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This procedure gave rise to what is now called the Abraham- Lorentz- Dime (ALD) self force in flat space [see Eq. 
( |l0| ) below for the analogous scalar particle case]. The concept of "radiative potential" was much later employed by 
Gal'tsov H for calculating the temporal rate of change of the energy and azimuthal angular momentum parameters 
for electrically charged particles orbiting a Kerr black hole. Though Gal'tsov analysis yielded correct results, it seems 
conceptually difficult, in general cases, to justify the use of such a non-causal approach. The problem becomes obvious 
when considering curved spacetime, where the self force exhibits nonlocal contributions: according to this approach, 
the force acting on a particle at a given moment turns out to be affected by the entire future evolution of the particle. 

A formal method for calculating the momentary self force in curved spacetime, which employs the purely retarded 
Green's function and is thus inherently causal, was developed long ago by DeWitt and Brehme Q (for the electromag- 
netic case). Dewitt and Brehme first carried out a local analysis of the retarded Green's function near the particle's 
worldline, based on the Hadamard expansion ^ (which is, basically, an expansion of the Green's function in powers 
of the geodesic distance between the source and evaluation points). Then, the particle's equation of motion was 
deduced by imposing local energy-momentum conservation on a thin world-tube surrounding the particle's worldline. 
Following the same approach (though using a different renormalization scheme), Mino, Sasaki, and Tanaka Q studied 
the gravitational self force by analyzing the local metric perturbations near a particle. They concluded that the reg- 
ularized gravitational self force in vacuum is due solely to a nonlocal tail contribution. It remained unclear, though, 
how to practically evaluate the formal expression derived for this tail (see, however, recent attempts to tackle this 
problem [|U)). 

Recently, Quinn and Wald Jll[] developed a different, axiomatic, approach to the regularization problem of the self 
force in curved spacetimes. Their approach relies on a "comparison axiom" , which allows the calculation of the self 
force by comparing the given problem with a (properly chosen) analogous problem in flat space (see for details). 
The implementation of this approach for both the electromagnetic and the gravitational cases |llj] yielded formal 
results in agreement with those obtained by DeWitt and Brehme (as corrected by Hobbs the main result by 
DeWitt and Brehme contained a trivial error) and by Mino et a/., with the advantage of involving much simpler 
calculations. More recently, the same approach was applied by Quinn [|13l for a scalar particle [the main results of 
this work, also quoted in Ref. 14 1, are summarized below; see Eqs. (P) - ([131)1 ■ Again, however, despite the availability 
of a simple formal framework for obtaining equations of motion for a test particle in curved spacetime, the practical 
implementation of the formalism in actual calculations (particularly, the evaluation of the nonlocal tail contribution) 
remained a challenging task. 

So far, the study of the self force effect in concrete situations have been restricted to very few simple cases. DeWitt 
and DeWitt |l5|] employed the above mentioned formalism by DeWitt and Brehme to study the self force correction to 
the geodesic equation for an electrically charged particle freely falling in a static weak gravitational field, in the limit 
of small velocity. They concluded that a repulsive force of cx q 2 r~ 3 magnitude (where r is the Schwarzschild radial 
coordinate) would be experienced by such a particle, in addition to the usual attractive inverse-square force. Later, 
Smith and Will [16| (and, independently, Frolov and Zel'nikov |l7|]) were able to derive an exact analytic expression 
for the 0(q 2 ) self force acting on an electrically charged particle held static in the Schwarzschild exterior. They 
found a repulsive self force of exact magnitude Mq 2 r~ 3 [measured by a momentarily static freely falling observer, and 
expressed in relativistic (geometrized) units] . This result was later extended by Frolov and Zel'nikov pq] to scalar and 
electrically-charged particles held static outside a charged, Reissner-Nordstrom type, black hole. It was concluded 
that, unlike in the electromagnetic case, no self force is experienced by a static scalar particle. This last result has 
been reproduced very recently by Wiseman |3) in a thorough analysis of the self force acting on a static scalar particle 
in Schwarzschild spacetime. 

The exact calculation in the static electric charge case was made possible owing to the existence of an exact analytic 
solution, discovered long ago by Copson [|l9| (and later corrected by Linet |2(|), for the electrostatic potential of a 
static charge in the Schwarzschild geometry. (The analogous closed-form solution for a scalar particle was constructed 
by Wiseman in S.) In more general cases one cannot benefit from the existence of such exact solutions. The usual 
approach for treating this problem, of obtaining solutions to the field equation in black hole backgrounds, is through 
the Fourier-multipole decomposition of the field. In the context of the radiation reaction problem this approach seems 
to offer two obvious advantages: First, it allows, in the usual manner, reduction of the field equation [originally a 
partial differential equation (DE) in 1+3 dimensions] to an ordinary DE, thus making it accessible to simple numerical 
treatment. Second, each individual mode of the field turns out to be continuous (and the corresponding self force 
to be bounded) even at the particle's location. Having this in mind, Ori previously proposed (21] that a practical 
calculation of the self force effect may be carried out by first evaluating the effect of each Fourier-multipole mode of 
a particle's self retarded field on its radiative evolution (through the local self force experienced by the particle), and 
then summing over all modes. 

In Ref. El}] the above sum-over-modes approach has been successfully applied for the calculation of the adiabatic, 
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orbit-integrated, evolution rate of the three constants of motion in Kerr spacetime, including the Carter constant.^ 
However, it appears that a naive application of this method for the calculation of the momentary self force would 
not be useful, in general. The reason is that, although each mode yields a finite contribution to the self force, the 
sum over all modes is found, in general, to diverge. This situation manifests itself already in the most simple case, 
that of a static scalar charge in flat space: in this basic example, the contribution of each multipole mode to the 
radial component of the self force is the same, —q 2 /{2r^) (where ro is the distance of the particle from the origin of 
coordinates, with respect to which the spherical harmonic functions are defined), with an obvious divergence of the 
sum over modes. This, however, does not mean that one has to abandon the whole sum-over-modes approach; one 
may still be able to benefit from its advantages, by introducing a suitable regularization procedure into the calculation, 
properly designed to overcome the above kind of divergence. 

In a previous paper j22| we introduced the basic elements of a method for the calculation of the self force in 
curved spacetime through regularization of the mode sum. The implementation of the proposed calculation method 
for a specific trajectory in a given spacetime consists of two stages. First, one solves (numerically, in most cases) 
the appropriate ordinary DE for each Fourier-harmonic mode l,m,ui of the retarded field, and evaluates the (finite) 
contribution of each of these modes to the self force. (Alternatively, one may numerically solve the 1+1 partial DE 
in the time domain, for each multipole mode l,m.) Then, the sum over all modes is made subject to a certain 
regularization procedure, which requires the knowledge of several regularization parameters. These parameters are 
derived analytically, for any given trajectory, through local perturbative analysis of the (retarded) Green's function. In 
Rcf. |2^| we outlined this regularization method as applied to a scalar particle moving on a Schwarzschild background, 
and presented final results (i.e., the values of all necessary regularization parameters) for the case of radial motion. The 
target of the present paper is three- fold: (i) providing a systematic presentation of the regularization scheme (including 
a discussion of some mathematical subtleties left untreated in (2^]); (ii) providing full details of the calculations involved 
in deriving the regularization parameters for radial trajectories; and (iii) extending the analysis to a wider class of 
static spherically-symmetric black hole spacetimes. 

This paper (as well as Ref. |2|]) is concerned with the analytic part of the regularization scheme; namely, it sets the 
mathematical foundation for the scheme, and demonstrates the calculation of the regularization parameters involved 
in its implementation (in the example of radial motion). As explained above, full calculation of the self force requires 
the supplementary numerical determination of the various modes' "bare" contributions to the self force. This was 
recently done for various trajectories of a scalar particle outside a Schwarzschild black hole: Burko first analyzed the 
case of a static particle |2i| and the one of a particle in circular motion |Q(see also ^5|). More recently, Barack 
and Burko applied the regularization scheme for studying radial trajectories in Schwarzschild pq |. These numerical 
works confirm the applicability of the regularization scheme, and provide support for the values of the analytically 
deduced regularization parameters. Of course, they also yield significant physical information. In the static scalar 
particle case, Burko recovered the familiar result, of a zero self force. Calculations of the self force on a scalar particle 
in circular and radial trajectories were carried out for the first time (see Refs. |^3|,^4],^6] for details). 

The current paper is organized as follows: In Sec. || we give some preliminary relations involving the self field, 
the Green's function, and the self force for a scalar particle. In Sec. [Ill] we decompose the Green's function into 
its spherical harmonic components, and discuss the applicability of this expansion. In Sec. ^ we decompose the 
(tail part of the) self force into its spherical harmonic contributions, discuss the need for regularization of the mode 
sum, and present the regularization scheme. The implementation of this scheme involves local analysis of the Green's 
function modes for large multipole numbers, which is carried out in Sec. |v|. The particular case of radial motion is 



then cons idered in Sees. VI and VII, where the regularization parameters for this case are explicitly calculated. In 



Sec. VIII we summarize, discuss possible extensions of the analysis, and briefly survey some related work. 



II. SELF FORCE ON A SCALAR CHARGE: PRELIMINARIES 

We consider a class of static spherically-symmetric (not necessarily vacuum) black hole geometries, having a line 
element of the formQ 

ds 2 = -f(r)dt 2 + f-\r)dr 2 +r 2 (d0 2 + sin 2 6 dip 2 ), (1) 



1 However, whereas the mode sum for the evolution rate of the energy and azimuthal angular momentum parameters was 
shown to converge (2l), it is not clear yet whether the corresponding mode sum for the Carter constant converges or not. 

2 Throughout this paper we use relativistic units (with G = c — 1), and metric signature — H — | — K 
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where t, r, 8, and ip are the Schwarzschild coordinates, and / is a function of r only, positive outside the event 
horizon. Important members of this class include the Schwarzschild solution, with / = 1 — 2M/r; the Reissner- 
Nordstrom solution, with / = 1 — 2M/r + Q/r 2 ; and the Schwarzschild-de Sitter solution, with / = 1 — 2M/r + Ar 2 /3. 
Here, M stands for the black hole's mass, Q represents its net electric charge, and A is the cosmological constant. 

We next consider a point-like particle of scalar charge q, with \q\ -C M, moving in a spacetime of the above type. 
Let x^ — x£(t) represent the particle's worldline (not necessarily a geodesic), with r being its proper time. The scalar 
particle exhibits a Klein-Gordon field $, satisfying 
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□ $ = $ ;a - Q = -±= y=gg a(> *, a ) )fi = -4mp{&), (2) 



where □ represents the covariant D'Alembertian operator, g is the metric determinant, and p{x^) is the scalar charge 
density, given by 

p(xn = q -=<5 4 (^-^(r))dr. (3) 
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The solution for the scalar field can be written as 

$(s") = J G^-x'^pix'^^d^x', (4) 

where G(x^\ x" 1 ) is the retarded Green's function, satisfying 

4-7T 

UG{x»;x m ) = -?=8 A {x> i -x"*), (5) 
V~9 

and subject to the causality condition, G = whenever x M lies outside the future light cone of x" 1 '. Combining Eqs. 
(||) and (^) we obtain for the scalar field 

/oo 
G[x»;x%(r)]dT. (6) 
-CO 

The "scalar force" experienced by the particle due to its own field shall be taken, following E3] , to be 

F a = g$ ;Q = q$, a , (7) 

evaluated at the particle's location. We comment that the so-defined force is not perpendicular to the four-velocity of 
the particle, u a = dXp(r)/dT, resulting in that the mass parameter of the particle is not conserved along the worldline. 
Indeed, the force on a scalar particle can be calculated otherwise (as in |6[, e.g.), such as to make the mass parameter 
conserved: = q($ t a + u a u^^^). Although we shall adopt here the simpler definition, Eq. (^), the results of our 
analysis could then easily be applied for the force F^ as well. (Given all vectorial components of F a , one can easily 
construct both the force component perpendicular to the worldline, and the component tangent to it.) With Eq. (g) 
we now have for the self force acting on the particle at a point Xq = x!^(t — 0) along its worldline, 

/oo 
G, a [x»;x%(T)]dT, (8) 
-oo 

where the gradient (taken with respect to x^) is to be evaluated at x M = Xq. 

The "bare" self force given in Eq. (||) needs to be regularized to avoid divergences associated with the behavior 
of the scalar field at the very location of the particle. For that goal, Quinn JlJ,[l4| applied the "comparison axiom" 
approach by Quinn and Wald for the scalar particle case. The total self force acting on the scalar particle was 
found to be composed of three parts: 

^p(total) _ ^.(ALD) _|_ ^(Ricci) _|_ ^i(tail)^ ^ 

The first term here is a local ALD-like term, reading 

i* ALD > = ^ 2 (a Q -aV), (10) 
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where a a is the four-acceleration of the particle, a 2 = apa^ , and an overdot represents covariant differentiation with 
respect to the particle's proper time r. The second term in Eq. (^) is related to the local Ricci curvature at the 
particle location, ft is given by 

^i Ricci) = \q 2 {R a pu^ + u a R M u p u-< - Ru a /2), (11) 

where R a p is the Ricci tensor and R is the curvature scalar. The third term in the expression for the total self force 
represents the non-local "tail" contribution. It may be expressed as 

F (taU) _ Um ( 12 ) 

where 

F^^q 2 ( ' G, a [x»;x$(T)]dT. (13) 

J — oo 

As we mentioned above, the occurrence of a tail term — a prominent feature of the self force in curved spacetimes — 
is due to the Green's function having its support also inside the source's future light cone. From the physical point of 
view, this is associated with the fact that waves are scattered off spacetime curvature while propagating on a curved 
background. 

The task of implementing the formal expression, Eq. (^|), in practical calculations of the self force is a challenging 
one. The difficulty stems, of course, from the need to evaluate the tail part, which requires knowledge of the Green's 
function everywhere along the particle's past worldline. Below we therefore focus on the tail term contribution to the 
self force, presenting a practical method for its calculation. 

III. MULTIPOLE DECOMPOSITION OF THE GREEN'S FUNCTION 

The regularization scheme to be introduced below is based on evaluating the contribution to the (tail part of the) 
self force due to each multipole mode of the (retarded) Green's function. To that end we first consider the multipolc 
decomposition of the Green's function. 

To begin, one may be tempted to decompose G into its multipole modes G l in the usual manner, as G — X)z°^o G l 
(where G l represents the quantity resulting from summing over azimuthal numbers to). Although this may look as 
standard procedure, caution is necessary here: in general, such a decomposition turns out ill-defined, as the sum over 
/ is found to diverge. This can be illustrated already in flat space. In this case, the modes G l admit a closed-form 
expression, which, for evaluation point lying inside the future light cone of the source point x , is given by 

c t (^,/) = (i±m^^il^. (14) 

Here, Pi is the Legendre polynomial, a = | Ut — t') 2 — (r — r') 2 ] , and 

cosx = cos cos 6' + sin0sin6''cos(i J 9 — ip'). (15) 



[Eq. ( |14|) can be verified by direct substitution, using Eqs. (|27D, (|28|), and (|30|), to be given below.] Consider, for 
example, the case x = 0, with P;(cosx) = 1 for all I, corresponding to both the source and evaluation points lying 
in the same radial direction. At large / values, the Legendre polynomial PiUt) admits the as ympt otic form cx Z _1 / 2 x 
oscillations with respect to I [the exact asymptotic form is given in App. [A| below; see Eq. ( A13[ )]. Thus, at large I 
one finds Gp lat cx I 1 / 2 * oscillations, implying that the infinite sum over all modes fails to converge]^] Mathematically 
speaking, this failure of the naive multipole decomposition may be associated with the fact that the Green's function 
of our problem exhibits a strong irregularity along the intersection of the future light cone of the source with the 
sphere of constant r and t. In Appendix |A| we discuss this irregularity in more detail, referring to the analysis by 
DeWitt and Brehme M . We then suggest a way to overcome the difficulty caused by the presence of the irregularity, 



3 In the more general case, with \ ^ 0, there is also aal x oscillations factor coming at large I from P;(cosx), yielding 
the asymptotic form Gp lat cx const x oscillations. Hence, clearly, the sum over modes G l diverges in the general case as well. 
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and explain how a well defined mode decomposition can still be accomplished. Although the detailed discussion of 
this issue is left to the Appendix, we outline here the basic argument, and present some definitions and notation 
needed in the sequel. 

In Appendix |A| we construct a "modified" Green's function G mo d = G — SG, where the function SG is chosen such 
that G mo d has the following properties: (i) it is a continuous function of and (p across the sphere of constant r 
and t; and (ii) it yields the same self force, through Eq. ([l3|), as the original Green's function G (this is guaranteed 
by taking the function SG to have no support inside the future light cone of the source). It is then argued that 
the modified Green's function admits an (absolutely) convergent multipole expansion, G mo d = J2i(G l — SG 1 ). Next, 

we define the new operation limj^oo ("tilde- limit") of a series of numbers Ai, as the standard limit lim^oo (when 

existing and finite) of the series Ai — B^ 1 ' — B^ 



(k) i / 1 . 

— -B ; , where the B t , s are any finite number of terms having 
the form B^ = cijl bj cos(ajl + /3j), with a,j, bj, ctj, and j3j being some Z-independent real numbers, and with none 
of the numbers aj vanishing. Namely, if there exist k quantities B\ of the above form, such that subtracting them 
from the original series A\ would yield a well-defined finite limit as I — > oo, then we define 



>0) 



linij. 



- 2 Ai = lim 

/ — >oc 



U) 



(16) 



We also define the "tilde-sum" of a series A\ by 

oo I 



(17) 



1=0 



1=0 



where ^ is the standard summation operation. It can be easily verified (see App. [a]) that when the "tilde-limit" 
(or the "tilde-sum") of a series exists, then it is unique. In particular, if a standard infinite sum Ai converges, 
then one may replace it with a "tilde-sum" operation. Thus, we may replace the convergent standard sum G mo d = 
J2i=o(G l ~ SG 1 ) with a tilde-summation. In Appendix [a| we show that J2i=o(^G l )=0. Consequently, we conclude 



G 



mod 



= E G - 



(18) 



1=0 



We emphasize once more that the modified Green's function G mo d can serve instead of the original function G for the 
calculation of self force, as both functions yield the same force. Eq. (|l8|) thus implies that the calculation of the self 
force can be carried out through analysis of the original Green's function's modes G l , by applying the tilde-summation 
instead of the (ill-defined) standard summation. (A more thorough discussion of the arguments leading to Eq. dlg| ) 
will be given in Appendix |a|.) 

Let us now turn to study the form of the multipole modes G l in greater detail. These modes can be written more 
explicitly as 



(19) 



=-! 



where Y lm {9, (p) are the standard spherical harmonic functions on the sphere of constant r and t. Substituting Eq. ( |l9| ) 
and the relation 6(9 — 0')5(<p— <p')/ s'mO = J2 l Yi m (0, ip)Y l * m (9' ', tp'), (where an asterisk denotes complex conjugation) 
in Eq. (p|), we obtain from the orthogonality of the spherical harmonics, 



r 2 r\r)g l ™ - [r 2 f(r)g 1 ™] _ p + 1(1 + l)g lm = 4n5(t t')5(r - r')YC m (0', 

In terms of the m-independent variable g (t, r;t' ,r'), defined through 

g lm = 2ng l Y l * m (9', l p')/(rr'), 

Eq. ( pp| ) becomes 

g% - + W l (r)~g l - 2f(r)S(t - t')5(r - r% 



(20) 

(21) 
(22) 



G 



where the radial coordinate r*(r) admits dr*/dr = f 1 (r), and the effective potential V (r) is given by 

v w =i m (>i+v + m\ (23) 



4 \ 7* V 

(with a prime denoting d/dr). 

To account for the causality condition, it is convenient to introduce the (Eddington-Finkelstein-like) null coordinates 

v = t + r* and u = t — r, . (24) 

The relation 8(t — t')5(r — r') = 2f~ 1 (r')5(v — v')5(u — u') can then be used to write Eq. (|22| ) in the simple form 

~g\ vu + V l {r)~g l = 8{v-v')8(u-u'). (25) 

We now impose causality by writing 

g l = g l (v, u] v', u')Q(v - v')Q(u - u'), (26) 
where Q is the standard step function. The "reduced Green's function" g (v, u; v' , u') obeys the homogeneous equation 

9 l ,vu+V\r)g l = (27) 



for all u > u' and v > v' . Substituting Eq. ( p6| ) into Eq. (25) and examining the behavior along the null rays v — v' 
and u = v! ', one finds that g l must admit 

g\v = v')=g l {u = u') = l. (28) 

For any fixed source point v',vf, the homogeneous equation (|27j), supplemented by the initial conditions (|2^), consti- 
tutes a characteristic initial- value problem for the function g^my where at u > v! and v > v' . 

Finally, to express G l in terms of the reduced Green's function g l , we substitute Eq. (|l]) [with Eq. (p6|)] into Eq. 
( |l9| ) . In the resulting expression we can explicitly sum over m by making use of the relation |2^| 

i 

J2 Y lm (9 1 <p)Y£ n {0',<p') = (An)- 1 (21 + l)P z (cos X ), (29) 

m— — I 

where cosx is the quantity given in Eq. (|l5|). We then find for the /-mode of the Green's function, 

G l = LP l (cosx) 9 ^ V,U 'f ,U ' ) Q(v - v') 0(u - u'), (30) 



where we have set 



L = 1 + 1/2. (31) 



IV. MODE SUM REGULARIZATION SCHEME 
A. The need for a mode sum regularization 



Following the discussion of the preceding section, we now replace G in Eq. (|13J) by G mo d, and then substitute for 
G mod from Eq. @. We finoQ 



4 It is assumed here that both the differentiation and the integration involved in constructing _Fi E ' out of G can be per- 
formed term-by-term with respect to the tilde-summation. This assumption should be verified by more closely inspecting the 
convergence properties of the tilde-sum over G l in Eq. (|l^), which, however, would be beyond the scope of the current paper. 
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^i £) = E^ (e) ' ( 32 ) 

where F 1 ^ represents the contribution to associated with the /-mode Green's function: 

F^^q 2 I e G\ a [x^ X ^r)]dr. (33) 

J — OO 

For practical reasons which become clear below, let us now write 

= f*-6F$ c \ (34) 

in which[] 

FL=q 2 G l a [<;^(r)] dr and SF^ = q 2 / G l a [a#;a£(r)] dr. (35) 



Here, F l a is the /-mode of F a = q& ia — the quantity given in Eq. (|^), which is sourced by the entire worldline. This 
quantity can be obtained from the /-mode of the self field, which, in turn, can be calculated essentially with no 
difficulty (using numerical methods, in most cases |23| , p4] , p6| ) . The other quantity appearing in Eq. (j34|), SFa e \ is 
local in nature, and thus may be treated, in principle, by means of local analytic methods (as we, indeed, demonstrate 
in this paper). 

In terms of F l a and 5Fa e \ the tail part of the self force is calculated through 



^(tail) 



oo 

^E^"^)' ( 36 ) 



To carry out this calculation, one may tempt to first calculate the sum over F l a (which is e-independent), and then 

evaluate the local contribution lim £ ^ + J2i^Fa^ ■ However, here one comes across a problem: Although each of the 
modes F l a yields a finite contribution at the particle's location, in general the sum over all modes F l a diverges. As 
we mentioned in the Introduction, this can be demonstrated even in the simple case of a static scalar charge in flat 
space. To overcome this type of divergence, the introduction of a certain regularization procedure for the mode sum 
is required. Such a procedure is described (and later implemented) in what follows. 

B. The regularization scheme 

To regularize the modes sum, one seeks a (simple as possible) e-independent function h a , such that the series 
^2i(F^ — h l a ) would converge. Once such a function is found, Eq. (|36|) can be written as 



-A., (37) 

1=0 

where 

56 

D - = 1™ E(^a (e) -*L)- (38) 
e_> 1=0 

In principle, a regularization function h a can be constructed by investigating the asymptotic behavior of F l a as 

I — * oo. It is also possible, however, to derive h l a from the large-/ asymptotic behavior of 5Fa . The latter and F l a 
must have the same singular behavior at the tilde-limit / — > oo (for fixed e), as their difference yields a convergent 



Strictly speaking, the two quantities and 5Fa are not well defined without specifying the direction through which the 



gradient of G is calculated. This issue is discussed in length later in this section. 
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tilde-sum over I. Obviously, in order to determine h l a (and D a ) from SFa, one merely needs the asymptotic behavior 
of S Fa ^ in the immediate neighborhood of e = 0. This allows one to derive h l a (and D a ) using local analytic methods, 
as shall be demonstrated in the next section. 

First, however, it would be necessary to comment here about a certain indefiniteness involved in the above definitions 
of the quantities F l a and SF& . 



Discontinuity of F l a and 5 Fa 



Whereas the quantity F a of Eq. (^3|) is well defined, the values of the two quantities F l a and 5F a depend on 
how exactly one evaluates the gradient G l a at the particle's location. To make this point clear, consider first the 

r-components Fj. and 6F^ e \ These are calculated according to Eq. (|35|) from the r-derivative of G l , reading 



G „, — 



£-Pz (cos %) 



{ [g] r - 9 l / r ] e ( v - w ')©(« - u') + / _ V [S(v - v')Q(u -u')-Q(v- v')5(u - u')] } 



(39) 



To calculate Fj and 5F r one needs to evaluate this derivative at the self force's evaluation point, = x$, with a 
source point x /M = x£(t). Now, if the derivative at Xq is calculated from [namely, by taking the limit r — 



G'(r)-G'(r ) 



], then the 8(u — u') term in Eq. 



-J of 

J39|) will have a nonvanishing contribution to Fj and to SFr [through 
v') term will have no contribution — see figure || On the other hand, 



the integrals in Eq. fl3q)1, whereas the S(v 

if the derivative is taken from r~ , it will be the 5(v — v') term to contribute, and the S(u — u') term to have no 
contribution. One can easily verify (as we explicitly do in the following section) that these two different 5 terms 

yield different contributions to the integrals in Eq. (Esh. Thus, although each of the quantities Fj and 5Fr has well 



defined values when calculated from either the limit 



or the limit 



'0 > 



these two 



nded values do not 



coincide. [Note that the quantity F r defined in Eq. (|33j) does not exhibit this kind of discontinuity, as for any finite 
e neither of the two 5 terms contribute to this quantity] 




Contribution 
from 5[v-Vp(i)] 



r ol r 
Contribution j 

from 5[u-J(p(T)] 



These quantities are calculated from Eq. 




by integrating over r the Green's function's 



FIG. 1. Discontinuity of F* and 5F" 
r-derivative at the particle's location, r = ro- If this r-derivative is calculated from the limit r — » rj (right figure), then the term in 
Eq. involving 6(u — u p ) will contribute to the r integration, whereas the 6(v — v p ) term will have no contribution. The situation is 

(left figure): then, contribution will come only from the S(v — v p ) term. 



reversed if the r-derivative is calculated from the limit r 
Since these two contributions are different, one finds that each of the quantities F}. and SF^S^ 
which are finite in magnitude). 



has two different one-sided values (both of 



One can similarly show that the i-components Fj: and 8F t also exhibit this kind of discontinuity through the 
particle's location (see the explicit calculation carried out in section [v] below). On the other hand, the 8 and tp 
components are obviously continuous through the particle's location, as G l depends on the angular coordinates only 
through the regular function P;(cosx)- 

For the sake of definiteness, we shall denote by F l + and SFr {e)+ the one-sided values arising from the r — > limit, 
and by F*~ and 5Fr~ the ones arising from the r — > limit. In addition, the symbols F^ and SF^^ will stand 
for the values derived from the limit t — > if dr/dt > at the force's evaluation point, or from the limit t — > t^ if 
dr/dt < there. (In case dr/dr = at the force's evaluation point, the two one-sided values of the ^-component turn 
out to coincide, as we obtain below.) With this notation we find, for a = r or t, 
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(and similarly for F^), where 



(40) 



Gf{x»rf») = [g\ r - 9 l /r T r 1 9 l 6{w ± - «/ ± )] , (41) 

and 

Gjft**;*"*) ee £J ^*> & + -«4)] ■ (42) 

Here, we have introduced the notation 

u; + = u and u;_ = v, (43) 

and, likewise, «/_ = u' and ml = v' . 



V. LOCAL ANALYSIS OF G l FOR LARGE 2 



The execution of the regularization procedure introduced above involves the construction of the quantities h a and 

D l a . As we pointed out earlier, this can be done by analyzing 5Fa E \ or, more accurately, the quantities SFa*^ given 
in Eq. (|4C|). For that goal, we must have sufficient information about the Green's function's /-mode G l , for large 
values of I, at the immediate vicinity of the self force evaluation point. In this section we use local analysis to obtain 
analytic approximation for G , up to the accuracy needed for the derivation of h l a and D a . 



A. perturbation analysis 



In Sec. ID above we have reduced the problem of calculating G l to that of solving a (l + l)-dimensional homogeneous 
partial DE for the function g l , Eq. (|2~7j), with the characteristic initial data specified in Eq. (pq). Given the function 
g l , the "four-dimensional" Green's function /-mode, G l , is then constructed from Eq. (|3C|). 

To explore the behavior of the function g l for small spacetime intervals and large I, we apply the following pertur- 



bation analysis. Let us separate the effective potential given in Eq. fl23| ) into two pieces, in the form 

V l (r) = L 2 V (r) + Vi(r), (44) 

where 

W = ^, and Vl (r) = ^[4rf(r)-1}. (45) 

Let us next expand V l (r) in a Taylor series in the small deviation r — ro about the particle's location r — ro- It is 
convenient to take the small expansion parameter to be r* — r*o, yielding 

V l (r) = V l {r ) + V\r )(n - r» ) + \v l {r ){r* - n ) 2 + • • • , (46) 

where an overbar denotes d/dr*. Let us also define 

A r = 2V o y 2 L(r* -r, ), (47) 

where Vbo = Vo( r o)- We shall refer to a variable of this kind, having the form L x (small spacetime deviation), as 
a "neutral" variable. Such "neutral" variables shall play an important role in our analysis, allowing one to properly 
take into account the delicate interplay between large I and small spacetime deviations. Expressing r* — r*o in terms 
of A r /L, substituting in the above Taylor expansion, and collecting terms of the same powers in L (with fixed A r ), 
Eq. ( f4q ) takes the form 

V\r) = V 00 [L 2 + L (/ x A r ) + (f 2 + f 3 A 2 r )] + 0(1/ L), (48) 
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where fx, fo, and are coefficients given by 



fx = \V^ /2 V = r 1/2 {rf - 2/), 
h = V^Vx = rf - 1/4, 

h = ~ V#o = y [(/') 2 // + /"] + 3(/ - rf). 



(Here, all quantities are evaluated at r = ro-) 

Denning now the dimensionless "neutral" coordinates 



y = Vqq ,L(v — v') and x = VqqL(u — u'), 



Eq. ( p7| ) becomes 



1+ fl^L + fl+j^t + (l/L s ) 



L 



L' 2 



g'=0. 



(49) 
(50) 

(51) 



(52) 



(53) 



We next expand g l in the form 



= J2 L k 9k(A r ,A r ,,z), 



k=0 



where the expansion coefficients gk are considered as being dependent of only the "neutral" variables A r , 



and 



A r , = 2V 1 ( ( 2 L{r[ - r, ), 



z = = (L/r )s. 



(54) 



(55) 



(56) 



Here, s is the geodesic distance, to leading order in r - r', between the Green's function evaluation and source points 
(when these two points have the same 6 and values): 



= {f(r )(v-v>)(u-u>)} 1/2 . 



(57) 



Substituting now the expansion ( |54| ) into Eq. ( p3| ) and comparing powers of i, we obtain a hierarchy of equations for 
the various functions gk, having the form 



9k,yx + 9k — Sk- 



(58) 



Here, the source Sk is determined for each fc > by the functions gw <k preceding gk in the hierarchy. In the analysis 
below we shall only need the terms with k — 0, 1, and 2. For these values of k, the source terms are given by 



S o = 0, 



(59a) 



Si — —fiA r g , 



(59b) 



S 2 = -fxA r gx-(f2 + f 3 A 2 r )g . 



(59c) 



Finally, to complete the formulation of a characteristic initial data problem for each of the functions gk, we supple- 
ment Eq. (|58| ) with the initial conditions 

9k(v = v') = g k {u = u) = S k0 , (60) 
which conform with the original initial conditions for g l , Eq. 
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B. Analytic solutions for k = 0, 1, and 2 

The solution to Eq. ( |58| ) for k = 0, subject to the initial conditions, Eq. (|60|), is given by 

.go = Jb(s), (61) 

where J„ are the Bessel functions of the first kind, of order n. 

To solve for gi, we first express the source Si explicitly as a function of y — x and z, using Eq. (|6l|) and the relation 
A r = y - x + A r > . We find 

Si = -fiJo(z)(y -x)- fiA r ,J (z). (62) 
Then, with the help of Table Q, we find the solution for k = 1 [satisfying Eq. (|60|)1 to read 

ffi = -|/izJi(z)(A r + A r O- (63) 
We now use the above solutions for go and 51 to express S 2 as 

S 2 = \f\zJi{z) [(y - x) 2 + 3A r ,(y - x) + 2A 2 r ,] - f 3 J (z) [(y - xf + 2A r ,(y - x) + A 2 r ,] - f 2 J (z). (64) 
With the help of Table |, we then construct the following solution for g 2 , satisfying Eq. ( |60| ) : 

.92 = ~zMz) [f 3 (Aj + A r A r , + A 2 r ,) + 3/2] + ^z 2 J 2 (z) [3f!(A r + A r ,) 2 - 8/3] + ^f?z 3 J 3 (z). (65) 



S{z;x- y) 


solution to g t y X + g — S(z; x — y) 


Mz) 

(y - x)Jq(z) 
(y - x) 2 J (z) 
zJi(z) 

(y - x)zJ 1 (z) 
{y-x) 2 zJi{z) 


zJi{z)/2 
{y-x)zJi(z)/4 

[z 2 J 2 (z) + 2(y~x) 2 zJ 1 (z)]/12 

z 2 J 2 (z)/4 

(y - x)z 2 J 2 (z)/6 

\z 3 J 3 (z) + 3(y-x) 2 z 2 J 2 (z)] /24 



TABLE I. Specific solutions to the inhomogeneous partial DE g iVX + g = S for various source functions of the form S(z; x — y), where 
z = These solutions can be verified by substitution. In this table, J n are the Bessel functions of the first kind, of order n. 



C. G t a expanded in powers of 1/L 

We are now in position to write the three leading-order terms in the 1/L expansion of th egradient G 1 ^[xq~, x%]. To 
that end we shall need, in view of Eqs. (^l|) and ([|2]), to evaluate the functions g^ derived above, along with their 
gradients gk,a, for x^ = Xq and x" 1 — x%. For the calculation of g^ a it is convenient to use the auxiliary relations 

d[z n J n {z)]/dz = z n J n _ 1 (z) (66) 

[for n = recall J-i(z) = — Ji(z)], along with 

dz/dr = — foL(y — x)/ z and dz/dt = ffoL(y + x)/z, (67) 

where 

/o = [rof^iro)}- 1 . (68) 



With the help of these relations wc derive from Eqs. (61 ), (63), and (pa) expressions for g^ r and gk : t (where k = 0, 1, 2) 



We then set x 11 = Xq and x = x£ in these expressions, and also in the expressions for the functions gk themselves 
(noticing the vanishing of A r ). All resulting expressions are then substituted in the formulas for G l fr and G l f, Eqs. 
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and (0). In these equations we also make the substitution S(w± — w' ± ) — 2LV^ 2 5{w±) where w± are 'neutral' 
variables defined by 

w + = 2x and w_ = 2y. (69) 
Finally, collecting common powers of L we obtain, for a = r or t, an expression of the form 

G[t K: - (^ 0)± i 2 + AWL + ( 2 > + •••), (70) 

where the various coefficients H a are functions of only the 'neutral' spacetime-interval variables^] A r /, z, w±, and 

A t = 2V 1 /2 L(i - t p ). (71) 

These coefficient functions are given by 

H( 0)± = fo[-A r >Ji(z)/zT Jo(z)S(w±)}, (72a) 

HP = —foft [zJ^z) + A 2 r ,J (z)] - J (z)/r , (72b) 

= Lf oAr , [7/ 2 z 2 J 2 (z) + 3zJ!(z)(/ 2 A 2 , - 8/ 3 ) - 16J„(z)(/ 3 A?, + 3/ 2 )] + -^A^zJ^z)/^ (72c) 

and 

#t (0)± = ffo [-AtJi(z)/z + J (z)S(w ± )} , (73a) 

H[ 1] =-\ffofiA r ,A t J (z), (73b) 

H[ 2) = ^//oA t [A 2 z 2 J 2 (z) + zJ!(z)(3/ 2 A 2 , - 8/ 3 ) - 16J (z)(/ 3 A 2 , + 3/ 2 )] (73c) 

(where the function / is to be evaluated at r = tq). In the above expressions for H, a) ' (2) and ff t (1) ' (2) we have omitted 
terms of the form oc z k J n (z)S(w±), with k being a positive integer, as such terms would yield vanishing contributions 
to 6 Fa e ' )± when integrated over r in Eq. (f4(i|). 

The quantities <5iv can now be constructed, in principle, for any given worldline, by inserting Eq. ( |70| ) into Eq. 
( pp| ) and carrying out the integration over r. In practice, to integrate over r, one should proceed as follows: Recalling 
that r is a small quantity (we have |r| < e), one first expands in powers of r all r-dependent quantities in the 
integrand of Eq. (|40[). (At that point, the details of the specific trajectory under consideration enter the calculation in 
an explicit manner; specifically, the power expansion coefficients turn out to depend on the values of u a , ii a , and u a 
at the force's evaluation point.) Then, since we are interested in extracting the large I (large L) behavior of SFa 6 ^, 
we introduce the "neutral" dimensionless proper time variable, defined by 

X=-(L/r )r, (74) 

and replace r by —ro(X/L). The integrand then takes the form of a power series in 1/L, with A-dependent coefficients. 
Transforming finally from integration over r to integration over A, one obtains an expression for 8F& in the form 
of a power series in 1/L, as desired. In the rest of this paper we carry out the above calculation in full detail (and 
derive h a and D a ) for the case of a purely radial trajectory. 



3 We hereafter use the symbols A r /, z, and w± to represent the values of these variables for x M = x^ and x' 
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VI. FORM OF THE REGULARIZATION FUNCTION h l : THE CASE OF RADIAL MOTION 



When considering radial trajectories (namely, ones along which d6 = d(f = 0) one has P;(cosx) = 1. Consequently, 
the Green's function given in Eq. (|30| ) becomes 9, ^-independent, resulting in the vanishing of both angular components 
of the self force (as should be expected, of course, by virtue of the background being spherically-symmetric). In the 
following we discuss the r and t components of the self-force. 

To carry out the integration in Eq. (|4^) we first expand each r-dependent quantity in the integrand in powers 
of 1/L, with A held fixed. The r-dependent quantities to be expanded are A r /, A t , 1/r', z, and the various Bessel 
functions appearing in Eqs. ( f72[) and (|7^). 

By expanding A r / in a Taylor series in r about r = rg, and transforming to the variable A, we obtain the expansion 



Aw = A r ,T + -A r /T 2 



p l/2 



-A r r, 
6 



1 

— L 

2 

-f.A + lh\ 2 (r /L) - \r\\ 3 (r /Lf 
2 o 



0(1/ Lf 



(75) 



where 



r* = f V, 

h = f- 2 (fr ~ f'r 2 ), 



■t\ = r 3 [(2/' 2 - /"/)r 3 - 3/'/rf + f 2 r] , (76) 
and where all quantities (except A) are evaluated at r = tq (t = 0). In a similar manner we obtain for A t , 

At = / 1/2 



i\-U\ 2 {r Q /L) + ^i\\r /L) 2 



and for 1/r', 



1 

r 7 



1 

ro 



l + r(X/L) + -(2r 2 -r r)(X/Lf) 



0{l/LY, 



0{l/L)\ 



Next, recalling z — (L/r )s, we obtain 



z = -s\ + bx 2 (r Q /L) - hx z {r /L) 2 + 0(1/L) 3 . 
2 b 



(77) 



(78) 



(79) 



To calculate the r-derivatives of s (which are understood here to be evaluated at t = 0), we make use of the 
normalization relation vu = 1//, and of the relations derived from it by successively differentiating its both sides with 
respect to r: vu + vu = (1/ f)'r, and vu + 2vu + vu = (l//)"f 2 + (l//)'f. Using these relations we find 



s(t = 0) 



-1. 



s(T = o) = -(/7/)f, 



s(r = 0) = 



16/ 2 



8/7-13/' ) r 2 +8/7r +-rfv* 



(80) 
(81) 

(82) 



(note that whereas r is non-positive throughout the integration domain, z and s are, by definition, non- negative) . 

Finally, we need to similarly expand the various Bessel functions appearing in the integrand of Eq. ([i"o|). Using Eqs. 
@ and @ we find for any n > 0, 

(83) 



Jn(«) = JnW + ^{r /L)sX 2 J' n {X) + (r /L) 2 ^-s'X'J^X) - -sX'J^X) 1 + 0(1/L 3 ), 

where a prime denotes d/dX. Using this general form together with Eq. (^6|), we obtain the following expansions, 
needed for our analysis: 
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JoO) = J (A) - -(ro/LYsX 2 ^) + 0(l/L 2 ) 
Ji(z) = MX) + ~(r /L)s(A 2 J (A) - AJ X (A)) 



-(r-o/i) 2 



;s 2 (A 3 J 2 (A) - A 4 Ji(A)) - -s (A 3 J (A) - A 2 Ji(A)) 



0(1/L 3 ). 



(84) 



(85) 



We now substitute the above expansi ons fo r A r / , A t , z, and the Bessel functions in Eqs. (|7^) (f73|). We also 
substitute for the delta functions in Eqs. (72a) and (73a) 



S(w±) 



S(X) 



S(X) 



>- = f^ 2 w T 6(X), 



(86) 



\dw±/dX\ / 1 /2 1 i, ± 

where the last equality is due to the normalization of the four-velocity. We thereby obtain expressions for the 
various functions Ha -0 ' 1 ' 2 ', each expanded in powers of 1/L up to order 0(L n ~ 2 ) (with A-dependent coefficients). 



Substitution of these expressions [and of the expansion for 1/r', Eq. (JTSj)] into Eq. ( |70D finally yields the desired 
expression for the Green's function's gradient, as a power series in 1/L (with A held fixed). We find (for a = r, t) 



G'± = H^L 2 + Hg>L + ff« + 0(1/L), 

(n) 

where the various coefficients H a are functions of A along the worldline, given by 



i 



[r*Ji(A)T^Jo(A)5(A)], 



(87) 



HP = -^3 [AJi(A) (/i// 1/2 - 4rr* + 4r Sr* + 2r r, 
+ A 2 J (A) (7 1/2 /ir 2 - 2r sn) + 4J (A) 



(89) 



and 



96r 3 



[48AJ (A) (2f - / 2 r* 



+4A 3 J (A) (3/rros// 172 + 6r 2 s 2 r, + 4r 2 sr* - 4// 3 f 3 + 6r 2 sr* - 6/ 1 / 2 f,r r,r. 



-12fr $f* + 6f 1/2 hrfl 



+8A 2 Ji(A) (3/if// 1/2 - 6r s - 3/ 3 r» - 12r 2 r„ + 6r TT* + 12r fsr* 

— 6r 2 s 2 r» — 4rg'sr* + 6r rf» — 6r 2 sf» — 2r 2 'r» + "if 1 ^ 2 fif^j 

+3A 4 J X (A) (4r 2 S 2 r* + // 2 r 3 - 4/ 1 /2 /l r oS r 2 ) 
+A 3 J 2 (A) (7/ 2 r„ - 12r 2 s 2 r,)] , 



(°)± = I. [MX)w T S(X) - Jr(A)i] 



H 



(i) 



4r 3 



|AJi(A) [-4ft + 2r (2St + t)] + A 2 J (A) (-2r St + / 1/2 /in<) } , 



(2) _ 



/ 

96rg 



{48AJ (A)/ 2 i 



(90) 
(91) 
(92) 



-4A 3 J (A) [r 2 (6s 2 i + 4si + 6si) - 3r (irsi + J 1/2 /i<r* + f 1/2 fiih) + 6f 1/2 firir* - 4// 3 ir 2 

+8A 2 Ji(A) (/ 3 f + 12f 2 t - 6r fi - 12r rsi + 6r 2 s 2 i + 4r 2 'si - 6r ri + 6r 2 si + 2r 2 'i' 

+3A 4 Ji(A) (-4r 2 s 2 < + 4/ 1 / 2 / 1 r stf, - // 2 tf 2 ) 
+A 3 J 2 (A)(-/ 2 i + 12r 2 s 2 i)}. 



(93) 
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Changing the integration variable in Eq. (Eoh from r to A, we now have 



5F^ = q 2 r 



Le/r _ 



LH^ + H^+H^/L + OiL- 2 ) 



dX. 



(94) 



The desired regularization function h l a is to be constructed such as to extract the large I singular behavior of SFa 
(while maintaining the simplest form possible) . By virtue of Eq. (S3) we take this function to be 



where 



h£ = LAf t + B a + CJL, 



(95) 



At = lim 



(96a) 



B a = lim 



I — »oo 



(5F l ^±~LA±)=q\ a 



(96b) 



C a = lim^ooi (SF*®* - LA± 



B, 



H^(X)dX, 



(96c) 



with J°°( )dX standing for lim^^oo f x ( )dX^\ That the second equality in each of Eqs. ( |96b| ) and (96c) is valid, and 
that the above choice of function h 1 ^ indeed satisfies the requirement that the tilde-sum J2h=o \SFa e ^ — h!j) would 
converge, will be shown in the next section, where we explicitly calculate the parameters A a , B a , and C' a . The reason 
for using the tilde-limit instead of the standard limit in the definitions of these parameters will also become clear 
then. As to the parameter D a , substituting now Eqs. ( |9^ ) and ( |95| ) into Eq. (|38|), one obtains^ 



D a 



-q 2 ro 



lim > 



1=0 J Le/r 



LH^ + H^+H^/L 



In conclusion, we find the tail part of the self force to be given by 



dX. 



(97) 



F (t al i) = J- (F*± - AlL -B a - CJL) D a , 



1=0 



(98) 



where, from the above construction of h l a , it follows that the sum over I converges at least as ~ l/l. The implementation 
of our regularization scheme thus amounts to analytically determining the regularization parameters A J , B a , C a , and 
D a , using Eqs. ( |96| ) and (|97|). For the calculation of the tail term, one may use either F l J (with A+) or Fj (with 
A~). Of course, one may also use any combination of these two one-sided quantities (e.g., their average). It should 
be emphasized here that the final result of the calculation, namely the tail term Fa^ (having a well defined value 
at the evaluation point), should be the same regardless of whether it is derived from one of the one-sided limits, or 
from the other, or, say, from their average. 

We finally point out that, although Eq. ( p8| ) has been developed here for radial motion, an expression of this form 
is also valid for any other trajectory ]29]] . The details of the specific trajectory under consideration would only affect 
the values of the various regularization parameters. 



7 Here we extend the definition of the tilde-limit, given in App. |a|, from discrete functions with index I to continuous functions 
of A. In an analogous manner, the tilde-limit of a function /(A) as A — > oo would be defined through the subtraction of a finite 
sum of functions of the form B^\X) = ajX bj cos(a_,A + (3j) (with ctj 7^ for all j). It is simple to verify that the tilde-limit of 
a function, when existing, is single- valued. 

8 We assume here that the contribution associated with the 0(L~ 2 ) term in Eq. ( [94] ) vanishes upon taking the limit e — » 0, as 
this contribution is of order O(e) (this becomes clear from the calculation of D a in the next section). However, a problem may 
occur if this term fails to yield a finite contribution when integrated over A. Here we do not further investigate the behavior of 
the 0(L~ 2 ) term, and just assume that the above potential problem is not realized. 
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VII. DERIVATION OF THE REGULARIZATION PARAMETERS FOR RADIAL MOTION 

To carry out the calculation of the regularization parameters in this section, we shall need the following integrals, 
the derivation of which will be described in Appendix |b| For k, n £ N we have 

-oc f (n + k- l)H/(n —k — 1)!!, < k < n, 

/ X k J n (X)dX = ^ (-l)( fe -")/ 2 (fc + n- l)!!(fe-n- 1)!!, even k-n > 0, (99) 
" [ 0, odd k - n > 

[in applying this formula for k = n, recall (—1)!! = 1]. If the difference k — n is a positive odd integer, then we also 
have 

/>56 

/ A fe J„(A)lnAdA = (-l)( fe -" +1 )/ 2 (/c + n-l)!!(fc-n-l)!!. (100) 
Jo 



A. Derivation of 

Substituting Eqs. (ff) and © into Eq. 
find, recalling Jo(0) = 1 and Jo(oo) = 0, 

Af = -j (r* =F w T ) = =F-j i, 

A±=/^(-i + ^)=±^f. (101) 
'o '0 

Note that the two one-sided values of A a are, in general, not the same. Consequently, as argued above, the function 
h a (and also SFa^) exhibits two different one-sided values. We note that the averaged value of the parameter A a , to 
be denoted by A a , is found to vanish: 

A a = ^(A+ + A-)=0. (102) 
This vanishing of A a seems to occur for all trajectories of a scalar particle, not only the radial ones considered here 



96a) and carrying out the tilde- integration [with the help of Eq. (B9j)] we 



B. Derivation of B a 



By the definition of B a in Eq. (36b) we have, after substituting for A^ from Eq. (|96a|) 



B a = -q 2 r a lim(^ c 



Le/r a 



+ q 2 r a 



H^(X)dX. 



(103) 



The first term here cancels out upon taking the tilde- limit I — > oo: At large / (and fixed e), each of the components 

Ha behaves as oc / _1 / 2 times oscillations with respect t o I [these components are linear combinations of Bessel 
functions, the asymptotic form of which is described in Eq. (|B6|)1. To leading order in I /I, this is also the form of the 



integral over Ha in Eq. (|l03j) (which is carried out over asymptotically large values of A). Hence, the expression in 
the squared brackets is found to diverge as oc l 1 ^ 2 times oscillations. When taking the tilde-limit, this divergent piece 
is removed, with the remaining part dying off at large I as oc Z -1 / 2 (ti mes oscillations). Therefore, no contribution 
arises from the first term in Eq. (103), and the second equality of Eq. ( |96b| ) is shown to be valid. 



To calculate the parameter B ai we now substitute for H a from Eqs. (|8S|) and (|92j) (for the r and t components, 
respectively). The calculation involves tilde-integrating over terms of the form oc AJi(A), oc A 2 Jo (A), and oc Jo (A). 
Reading the values of the these integrals from Eq. (|99|), and substituting for fi, r*, f*, and s [using Eqs. jis|), (|7q), 
and (pi])], we obtain 
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B r = -$ 1 r 1 {f + r f/2-r 2 +r f) (104) 

and 

B t = £[fr i+ir(r f'-f)]. (105) 

Recalling a r = r + f'/2 and a* = t + in the case of radial motion considered here, we may express this result in 
a more compact form as 

q 2 

B a = - (5£ + roa a - fu a ) . (106) 
Zr Q 



C. Derivation of C a 



By its definition in Eq. (96c), we have for the parameter C a , after substituting for A a and B a , 



Ln/ro 



(107) 



Again, there is a residual part left from the calculation of A Q and _B Q , which involves integration over asymptotically 
large values of A. This part can again be shown to vanish as the tilde- limit I — > oo is taken, resulting in that only the 
second integral in Eq. (107) survives. 

To calculate C a , we thus use the second equality of Eq. (96c), in which we substitute for from Eqs. ( po|) and 



( |93| ) (for the r and t components, respectively). One then has to evaluates the tilde- integral of a sum of various terms 
of the form oc A fc J„(A), all of which have k > n and odd k — n. According to Eq. (|99|), all such integrals vanish. Thus, 
we find 



C a = 0. 



(108) 



The vanishing of the parameter C a seems to be a universal feature of our scheme, regardless of the specific trajectory 
under consideration |29(| . As we also find below, this vanishing constitutes a necessary condition for the self-consistency 
of the whole regularization scheme. 



D. Derivation of D a 

To calculate the parameter D a we write Eq. (|97]) in the form 

D a =DW+D£) + D%\ 

where 



= -g 2 r 



hm V / L 1 - n H^ ± (X)dX. 



(109) 



(110) 



In calculating the above three pieces of D a , we shall transform from summation over I to integration over a continuous 
variable. For this transformation we will make use of the relation 



oo „ c 
J2(e/r )K(Le/r ) = / 

l_n JO 



K(x)-±(e/r fK"(x) 



dx + 0(e 3 ), 



(111) 



where K(x) is any (sufficiently regular) integrable real function, and x is an integration variable. Here, the O(e ) 
term on the right-hand side (RHS) is the standard ("Riemann type") integral, which is obtained, by definition, when 
the e — >• limit of the left-hand side is taken. We also indicated here the 0(e 2 ) correction to the integral, which we 
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shall have to take into account in the calculation below (it is straightforward to verify the form of this correction 
term using standard calculus). Obviously, Eq. (Ill) also holds for the tilde-sum YliZo' wnere on the RHS we use the 
tilde- integral J °°. 

Beginning with the calculation of , we write Eq. ( |l!0|) for n = as 



= -q 2 r Km 



(r / e ) 2 ^(e/r ) (Ze/r ) / H a °^(X)dX 

1=0 \ JLe/ro 



Comparing the form of the sum in this expression to the left-hand side of Eq. ( Ill ), we find 



= -q 2 r lim 



(r /efj Q °°dx (kW(x) - ^(e/r ) 2 [K^ (*)]") 



(112) 



(113) 



where 



kV>(x)=x H a °^(X)dX. 



(114) 



Note here how the contribution to due to the 0(e :i ) term appearing in Eq. (Ill) vanishes upon taking the limit 
e — > 0. We recall that each of the two components Hr and ifj , given explicitly in Eqs. (p8[) and (pl|), contains 



two terms: one proportional to <5(A) and the other to Ji(A). The oc (5(A) term has no contribution to K a , resulting 
in that both two one-sided values H a and H a ' yield the same function K a ^ (for that reason, no ± sign has bee n 
assigned to this quantity). There is an apparent danger of divergence coming from the 0(e -2 ) term in Eq. (113). 
Such a divergence is avoided, however, as we have 



K a °\x)dx 



'0 -'0 



dxx I Ji(X)dX = — 3 



x 2 Ji(x)dx = 0, 



(115) 



o Jo 



where in the second equality we integrated by parts with respect to x, and where the vanishing of the last integral is 
implied by Eq. (99). Thus, the 0(e -2 ) term in Eq. (112) vanishes, and the e — > limit in this equation turns out well 
defined. The remaining 0(e ) contribution reads 



^ 0) = ^?/ [K a °\x)]"dx = - 



24 



q 2 rg 
24 



H^ ± (x)+\xH a °^(a 



dx = 



q 2 r 

24 



H a °^(x)dx, (116) 



as the surface term vanishes. Substituting for Ha ± (x) and integrating using Eq. (|99j), we finally obtain 



M 0) 



q 2 u a 

' 24r 2 J(l 



Jx(X)dX 



q 2 u a 
' 24r 2 



We next turn to calculate Da' . Writing Eq. (Illfj) for n = 1 in the form 



L>W = -q 2 r lim 



°° poo 

(r /e)J2(e/r ) / H a ^(X)dX 

1=0 J Le/r 



(117) 



(118) 



and applying the summation formula (111), we obtain 

D a 1] = -q 2 r lim (r /e) / dx j H a 1 \X)dX 



OO pOO 



(119) 



with all 0(e 2 ) terms appearing in Eq. (Ill) vanishing in the limit e — ► 0. By integrating the last expression by parts 
with respect to x (noticing the vanishing of the surface term) one finds 



D a 1] = -q 2 r a lim 



(r /e) / xH a l \x)dx 



(120) 
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Here, the integrand contains only terms of the form oc x k J n (x), with k — n being positive odd integers [see Eqs. (p 
and (f)^)]. Hence, by virtue of Eq. (^9|), the integral vanishes, yielding 

= 0. 



Finally, from Eq. (11C) with n = 2 we obtain 



= -q 2 r lim 
e— >Q 



00 POO poo poo 

J2(z/ro)(Le/r )-i H™(\)d\ = -q 2 r Q (dx/x) H^(X)dX, 

2_q J Le/ro **0 Jx 



which, after integrating by parts, becomes 



D (2) 



-q 2 ro 



poo / poo ' 

J In x H {2) (x)dx - Urn ( In x J H (2) (XjdX 



(121) 



(122) 



(123) 



We notice here that the second integral on the RHS is just C a + 0{x) [up to a multiplicative constant; see Eq. ( 96c)]. 
The above-deduced vanishing of the parameter C a guarantees the definiteness of the x — ► limit in Eq. (123), and 
makes the second term on its RHS vanish. Note the way the vanishing of the parameter C a appears as a necessary 
self-consistency condition in our scheme: had we got C a 7^ 0, the parameter D a would have been indefinite, and the 
whole regularization scheme woul d have been rendered meaningless. 
As the x — > limit in Eq. (123) vanishes, we are left with 

= —q 2 ro / \nxH {2) {x)dx. 



(124) 
Jo 

With the explicit form of H a 2 \ given in Eqs. (^) and (|9^), the integrand in the last expression is found to consist of 
various terms of the form oc x k \nx J n (x), with k — n being positive odd integers. The integrals J °° of such terms can 
be read from the formula (10C), which we derive in Appendix [b|. Using this formula we obtain expressions for D {2) 

(2) 

and D t , which, after substituting for r*, r* 3 f*, s, s\ /i, and /3, read 



/(4/ - 3) + 8r /'(/ - r f) + 2rg(3//" - f) 



(125) 



and 



D 



(2) 



-\q 2 {fr + fHvu)- l -q 2 f>ri 

+ 2§ [3/ " 4/2 + WW " 8 W ~ ±f'r 2 o? ~ 2 //' V o ~ 8//V0] 



(126) 



We are now in position to write an expression for the "overall" parameter D a . We have D a — + D a z> , yielding 



,(2) 



D r 



1 



79 2 (/ r + ™«) + 



and 



2 • 
q r 

1 



2/(/ - 1) + 4r /'(/ - r r) + rg(3//" - / ,2 ) 



(127) 



A = -g9 2 (fr+f 2 ivu)~-q 2 f'ri 

+ 7?2 t 2 ^ 1 - /) + 2/-l(/') W - 4/'V 2 r 2 - 2/'rgf - ff'r 2 - 4//V ] 



(128) 



i/' and a* 



In the case of radial motion, the four-acceleration's components admit the explicit form a r = r 
t+f'f*t. Recalling also that in the spacetime class considered here the Ricci scalar reads R = — [/"+4/'/r+2(/— l)/r 2 ], 
one can show that the above two expressions for D r and D t can be put into the simple vectorial form 



D a = T^q 2 (a a - a 2 u a ) ~ ^q 2 Ru a . 



(129) 

Comparing now this result with Eqs. (]1C|) and (|llj), and recalling that in the case considered here (that of radial 
motion on static spherically-symmetric background) the first two terms in the expression for F a Ricci) [Eq. cancel 
out, we arrive at the remarkable conclusion that D a is exactly the "standard" local part of the self-force: 

D a = Fi ALD *> + Fi Ricci \ (130) 
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VIII. SUMMARY AND CONCLUDING REMARKS 



The total self force actin g on the radially moving scalar particle is obtained by substituting from Eq. ( |98| ) in 



Eq. (|9j). By virtue of Eq. (130), the contribution of D a to the tail term is then found to exactly cancel out the local 
term in the expression for the total self force. This, in addition to the vanishing of the parameter C a , leads to the 
simple result 

oc 

^totai) = J- (F l ± - AtL - B a ) (131) 



1=0 



(where, we recall, L = I + 1/2). An even simpler form is obtained when calculating p^° tais > usin g the averaged value 



of the modes F l a) obtained by averaging over their two one-sided values. Then, by virtue of Eq. (102), we find 



F (totai) =J2(F l a - B a ) , (132) 



1=0 



where F l a = ^(F^ + F l a ). Recall that the parameter B a , given in Eq. ( |106| ), is just the asymptotic value of the 
averaged l-mode F l a at the limit I — ► oo; namely, the total self force is obtained by simply subtracting from each (two 
one-sided averaged) mode its large I asymptotic value, and then summing over all modesj^j 

To summarize, in this paper we have developed a method for calculating the self force on a scalar particle in 
curved spacetime, through regularization of the multipole mode sum. The basic difficulty in applying the mode 
decomposition approach — the apparent divergence of the sum over modes — has been taken care of by the introduction 
of an appropriate regularization scheme, providing a practical prescription for calculating the self force. It should 
be emphasized that the proposed method does not involve any weak-field or slow- motion approximations, and thus 
allows effective calculations of the self force even for strong field orbits. 

The basic expression for the tail part of the self force is given in Eq. (|98|) . This expression was developed here for 
radial motion; however, the same general form applies for any trajectory |29(] , with the details of the orbit encoded only 
in the values of the regularization parameters (as well as, of course, in the form of the "bare" modes F^). To apply 
this general expression for a given trajectory requires knowledge of (i) the modes F l a , to be derived by supplementary 
(basically straightforward) numerical analysis, as done in Refs. |^,^J,^6|; and (ii) four regularization parameters for 
each spacetime component of the force. In this paper we have worked out the entire calculation of the regularization 
parameters for the case of radial motion. For any other trajectory, the derivation of these parameters can be carried 
out along the same lines, based on the explicit form of the Green's function's l-mode given in Eq. (|30|), with Eqs. 
(|6l|), (|63|). and (|65|). To that end, one first obtains an expression for the gradient of the Green's function's Z-mode, as 
in Eq. (|70|) [supplemented with Eqs. ( f72|) and ([T3])]. One next expands this gradient in powers of proper time r along 
the worldline about the force's evaluation point, and re-expresses the resulting expansion as an expansion in powers 
of 1/L, by holding tL fixed — as in Eq. ( p7\ ) [supplemented with Eqs. (^8|)-(P3|)1. One finally uses the values of the 
above expansion coefficients [denoted in this paper by Ha (tL)} to construct the regularization parameters through 
Eqs. @ and @. 

For a scalar particle movin g radially in a spacetime of the class considered in this paper, we found the total self 



force to be given by Eq. (132). This constitutes our main result for the radial motion case, together with the explicit 



values of the regularization parameters given in Eqs. (1021), (|106|), (108h, and (129). We have found that, in the radial 



motion case, the parameters A a and C a both vanish [the one-sided values of A a do not vanish; they are given in 



Eq. (101)], and the parameter D a is just the standard local part of the self force. The vanishing of C a , shown here 



explicitly, appeared as a necessary condition for the dcfiniteness of the whole scheme (had C a ^ 0, the parameter 



D a would have diverge; see the discussion in Sec. VII). This point serves to demonstrate the self-consistency of the 
regularization scheme. 

A question arises, whether the above results (the vanishing of A a and C a and the special value of D a ) represent 
generic features of the regularization scheme, or rather are special to radial trajectories. Preliminary investigation |29| ] 
suggests that, indeed, A a and C a vanish for all trajectories, at least in the Schwarzschild case. As to the parameter 



9 The simplicity of our main result, Eq. (132), may lead one to wonder whether there could be simple arguments leading 
directly to this result. Such arguments might perhaps rely on general properties of the Hadamard expansion. This should make 
an interesting subject for further investigation. 
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D a , this was shown so far to o bey Eq. ( J13Q ) at least in one more important example, that of a circular orbit ar ound 
a Schwarzschild black hole [t?9|j22] l . It might be conjectured (and be subject to further investigation) th at E q. (|l30| ) 
holds for any trajectory in any static spherically-symmetric background. In that case, the simple Eq. ( |132| ) for the 
total self force would be valid for all such trajectories. 

Under the above conjecture, we find that regularization of the total self force requires knowledge of just one 
parameter, B a , representing the asymptotic value of the modes as I — ► oo. Th e val ue of B a for any specific radial 



trajectory on any given spherically-symmetric spacetime, can be read from Eq. (106) (valid regardless of the above 



conjecture). For example, in the special case of a static particle we find B[ static) = and 



B (rtatfc) =_!_ (l + ^Jl ). (133) 



For radial geodesic motion we find 



^(geodesic) = _jf_- E ^ ^(geodesic) = ^ _ £ 2// ) > (m) 

2r 2r Q 

where E = —Ut is the energy parameter (which is a constant of motion in the absence of self force effect). For the value 
of B r in the case of uniform circular motion in Schwarzschild spacetime (the derivation of which will be presented 
elsewhere |2!| ) we refer the reader to Eq. (34) of Ref. p2| . 

The applicability of the regularization prescription described here was demonstrated recently in actual calculations 
of the self force for various scenarios. Burko first studied the cases of static |^3| and circular H] orbits in the 
Schwarzschild spacetime. For these stationary scenarios, the modes F l a were obtained by summing over the Fourier- 
multipole modes i 7 ^" 1 ", first derived by solving the appropriate ordinary field equations in the frequency domain. 
Later, Barack and Burko |26f| analyzed the case of radial motion in Schwarzschild. In this case, which is no longer 
stationary, numerical evolution of the appropriate partial DE in the time domain was applied to directly infer the 
modes F l a . In each of these studies, the overall force acting on the scalar particle was finally deduced by summing over 
all modes, using the above regularization scheme. In each of the cases analyzed, the vanishing of C a was demonstrated, 
and the analytically-derived expressions for A a and B a were verified. Aside from demonstrating the applicability of 
the regularization scheme and providing verification for the values of the regularization parameters, the above studies 
yielded valuable physical results, as mentioned in the Introduction. 

Of course, the analysis of the scalar self force merely serves as a toy model for more realistic cases. Generalization 
of the regularization scheme to the electromagnetic self force seems possible, based on the existing formalism [Q,ll|. 
Such a generalization is necessary, for example, for resolving the interesting question re-raised recently by Hubeny [ 30 1 , 
whether a nearly extreme electrically-charged black hole might be overcharged (and its event horizon by destroyed) by 
throwing in a charged particle: as pointed out by Hubeny, knowledge of the exact radiation reaction effect is crucial 
for obtaining a definite answer. More difficult to accomplish would be the important generalization of the scheme to 
the gravitational self force acting on a mass particle. 

Finally, it should be mentioned that a closely related approach was recently applied by Lousto |3l[] for studying the 
gravitational self force on a mass particle in Schwarzschild spacetime. This approach is also based on the multipole 
expansion, yet it employs a different regularization method for the mode sum (it is argued that the correct self force 
can be deduced by applying the zeta- function regularization technique). For the geodesic motion case studied by 
Lousto, this approach leads to an expression analogous to Eq. (|132|). 
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APPENDIX A: MULTIPOLE EXPANSION OF THE GREEN'S FUNCTION USING 

"TILDE-SUMMATION" 

DeWitt and Brehme wrote a general expression for the scalar Green's function in curved spacetime [see Eq. (2.21) 
of Ref. @], of the form 

Gix^yx" 1 ) = a(x";x^) 5(a) + b(x^;x'^) 6(a). (Al) 
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Here, a is plus or minus half the squared geodesic distance between the source point x^ and the evaluation point x , 
according to whether the geodesic connecting the points (along which the invariant distance is measured) is timelike 
(a > 0) or else (a < 0); a is a certain function having a well defined value at a = 0; and & is a function which may be 
written as a Taylor expansion in a about er = 0. [This expansion was shown by Hadamard (see pp. 96- 98 in Ref. (^]) 



to converge uniformly at least inside the region where a is single valued.] While the first term in Eq. ( Al) is associated 
with the familiar delta function exhibited already in flat spacetime, the second term represents a curvature-induced 
tail, which "fills" the light cone (defined by a = 0). Note that the Green's function is strongly irregular along the 
light cone of the source point x . 

Now, the question to consider is whether the above Green's function may be expanded in terms of the standard 
spherical harmonic functions Y lm (9, ip) on a sphere of constant r, t. Standard theorem (see, e.g., Ref. |32|, p. 513) 



states that a sufficient condition for (absolute and uniform) convergence of the spherical harmonic expansion is the 
expanded function being G 2 on the sphere. This condition is not satisfied here, as the Green's function diverges along 
the curve generated by the intersection of the future light cone of x ni and the sphere of constant r, t. Therefore, it is 
not guaranteed, in advance, that such an expansion could be naively applied. Indeed, already in flat spacetime the 
attempt to apply the multipole expansi on to the Green's function turns out to yield a divergent sum — see Eq. (|l4| ) 



and the discussion proceeding it in Sec. Ill 



Let us introduce the "modified" Green's function (which is not a "Green's function" for the scalar field anymore), 
defined by 

G mod = G-SG, (A2) 

where 

SG = a S(a) -b @{-a), (A3) 

with ao = a{<r — 0) and bo = b(a = 0). The function G mo d has two essential features: (i) its support inside the light 
cone is identical to that of the Green's function G (as SG has no support there), and (ii) it is continuous throughout 
any sphere of constant r, t. The first feature implies that we can use G mo d instead of G in calculating the self force: 
We may re-write Eq. (|l3|) as 

F« ee g 2 [ ' {G mod [a£; a£(r)] } Q dr, (A4) 

J — oc 

as the SG term in G mo d contributes nothing to the integral along the particle's worldline. 

The second of the above features of G mo d, its continuity, may imply that the multipole expansion could now 
be applied to it. Strictly speaking, G mo d does not satisfy the above sufficient condition for absolute and uniform 
convergence of the mode sum (i.e., being G 2 on the sphere); yet, we shall assume here that G mo d, being continuous, is 
already regular enough to admit a convergent mode sum. The results of our analysis turn out to be consistent with this 
assumption, as the mode-sum of G considered (based on the function G mo d through the use of "tilde summation" — see 
below) is found to be (absolutely) convergent. The validity of this assumption can also be d emonstrated in the flat 
space case: the oc bo tail term vanishes in this case, and by expanding the S(a) term in Eq . ( [A3] ) in spherical harmonics 
one can easily verify that the /-mode of SG is exactly the l-mode of G, given in Eq. (|1J). The l-mode of G mo d then 
vanishes, and the mode sum converges. Although trivial, this flat-space example may serve to demonstrate how the 
subtraction of 5G from G already removes the divergent piece from the l-mode, making the mode sum converge .p| 
We thus expand G mo d as 

oc oc 

G mod = Cd = E (° l - 6Gl ) > ( A5 ) 

1=0 1=0 

where G mo J and SG 1 are the spherical harmonic modes of G mo d and <5G, respectively (obtained by summing over all 
azimuthal numbers m). 

To proceed, let us now define the new operation lim^oo ("tilde-limit") as follows: Consider a series of numbers A\ 
(with / = 0, 1, . . . , oo). Let be any expression of the form 



10 One may similarly construct a more sophisticated function SG, designed to yield a C 2 modified function G mo d [by canceling 
also the O(a) and 0(a 2 ) terms in the Taylor expansion of b] being sufficiently regular to assure uniform and absolute convergence 
of the multipole expansion, by standard mathematical theorem. Such an improved construction will not be examined here. 
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B 



00 



ajl bj cos(ajl + (3j), 



(A6) 



where cij, bj, otj, and @j are some ^-independent real numbers, with a,j 7^ for all j. If there exists a finite number 

k of expressions B^ of this form (with j — 1,2,..., k), such that subtracting their sum from the original series A\ 

would yield a well-defined finite limit as I — > 00, then we define the "tilde-limit" limj^ooA; as in Eq. (|l6|). 

One may easily be convinced that the tilde- limit is single- valued (when existing). For, suppose that for a given 



series Ai there were two different sets of quantities B*f' , one (denoted by B\ 3> ) yielding lim;^ooA; = c, and the other 
(denoted by b\^) yielding limi_,oo Ai = c ^ c. Then, for the difference between the two limits one would have found 



Uh 



lim/_ 



>00 



00 



c-c^ 0. This, however, is impossible, as the (standard) limit I — > 00 of any quantity 



of the type £> ; is either diverging or zero, and so is the limit of any finite sum of such quantities. Hence, we must 
have c = c, and the tilde-limit is single valued. In particular, we find that if there exists a finite standard limit 



lim; 



Ai, then lim, 



Ai = lim, 



-4, 



We can now also define the "tilde-sum" of a series Ai, as in Eq. (|17|). Again, if the "tilde-sum" of a series exists, 
then it is unique. Also, if the standard sum X^z^oo converges, th en w e may replace it with a "tilde-sum" operation. 
In particular, we may replace the convergent standard sum of Eq. (|A5[) with a tilde-summation: 



G m od — ^ (G l — SG 1 



Below we show that 



1=0 



(A7) 



J2 SG l = 



1=0 



(A8) 



(for any evaluation point Xq lying inside the future light cone of the source point x%). As a consequence, the tilde-sum 
of G l is found to be finite and equal to G mo d — as indicated in Eq. (|l8|) of Sec. |Fn. 
Combining Eqs. (A4) and 



we conclude that the self-force can be calculated by analyzing the modes G l of the 
original Green's function, provided that in order to sum over all modes one applies the i^de-summation instead of 
the standard summation. The validity of this statement crucially depends on the vanishing of the tilde-sum over 8G l 
[Eq. (|A8|)], which we now prove. 



Proof of Eq. (A.8) 



To calculate the Z-modes of SG, it is convenient to use a spherical coordinate system in which the source point a;^ 
lies on the polar axis (i.e., 9 p = 0). In this coordinate system, contributions to SG would come only from the m = 
modes: 



5G l = V Y lm (6^) [ dip 1 [ d(coae f )5G(a , )Y^{e',< f /) 
, , Jo J-i 



m=—l 
pi 



d(cos0') [a S(a') - b Q(-a')] Pi (cos 9) Pi (cos 9') 



(A9) 



where the integration is carried out over a sphere spanned by 9' ' , ip' (containing the evaluation point ir M ), L = (Z + l/2), 
and a' is half the squared geodesic distance between the source point and the integration point 9',ip'. Now, the future 
light cone of the source poi nt (along which a 1 = 0) intersects the integration sphere along a circle 9' = const = 9q. 
The S(a') term in Eq. (A9) contributes to the integration only along this circle, while the 8(— a) term contributes 
only across the part of the sphere outside it. We thus find 



SG 1 = LPi (cos ( 



a Pi{cos9o) - b 



P/(cos 9')d(cos 9') 



SG 



(5) 



SG 



(e)> 



(A10) 



where do = 0,0/ \da' /d(cos9')\ g , =go is a constant, and where the symbols 5G l s and SGq represent the two terms 
proportional to do and bo, respectively. 
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To calculate the tilde-sum over SG l we make use of the finite-sum identity [referred to as "Christoffers first sum- 
mation formula"— see, e.g., Eq. 3.8(20) of Ref. |3] 



' 1 
(x - y) LP v {x)P v {y) = - (I + 1) [P t+1 (x)Pi(y) - Pi +1 (y)Pi(x)} 



(All) 



Z'=0 



valid for \x\ < 1 and \y\ < 1. Applying this formula to Eq. (A10), we obtain 



E K 



i'=0 



- a lim ; _ 



(-5) 



SG 



(e; 



P ;+1 (x)P i (x )-P i+ i(x )P ; (x) 



X — Xq 



1 



6 lim ( 



a+i) 



P i+1 (x)P(x')-.P;+i(z')^(z) 



(A12) 



where x and xq stand for cos 9 and cos#o- (Recall that for any evaluation point x M lying inside the future light cone 
of the source point xg, we have cos 9 > cos^o > cos 9'; hence, the denominators appearing in the last expression are 
strictly positive.) 

Let us consider first the tilde-sum of SG, g s . For values of 9 satisfying e < 9 < n — e (where e > 0) we have the 
large-Z asymptotic form 



Pi (cos 6) oc V2(Z7rsin(9r 1/2 cos(Z0 + tt/4) + 0(r 3/2 ) 



(A13) 



[see, e.g., Eq. 3.9(2) of Ref. p3]| ]. Using this asymptotic form with Eq. ( A12| ), we may easily write Y2i'=o ^ G \s) as a sum 

of a few terms of the form a,j cos{oijl+ (3j) + 0{l~ x ) (in the case < 9 < it), or of the form ajl 1 / 2 cos(ajl+ f3j) + 0(l^ 1 ^ 2 ) 
(in the case 9 — or tt), where a,-, aj ^ 0, and /3L- are certain functions of and # (independent of £). Such terms all 
vanish at the tilde-limit I — ► oo. Thus, we clearly have 



(A14) 



We next turn to calculate the tilde-sum of 6Gl e y Integrating by parts in Eq. ( A12| ), and using 



Pi (x) dx = 



Pi +l {x) - Pi^jx) 
21 + 1 



[see Eq. 7.111 of Ref. p7j, together with Eq. 8.733-4 therein], we obtain 



(A15) 



(e) 



l'=0 



x — x' 



-1 



X P^xjx^-P^xjx') ■ / 
[x — X ) z 



21+3 



Pi +2 (x')-Pi(x') 



rx 
J-l 



{x — x') 2 



dx' 



(A16) 



The surface terms here vanish at the lower boundary, as p + i(— 1) = P;_i(— 1) and P/+2(— 1) = p(— 1) [recalling 



oc 



P/(— 1) = (— 1)']. The contribution from the upper boundary dies off at large I [by virtue of Eq. ( A13[ )] as 
(times oscillations) for 9 ^ 0, or as oc l^ 1 / 2 (times oscillations) for 9 = 0. In both cases, this contribution thus 
vanishes at the standard limit I — > oo, and hence also at the tilde-limit. Consider next the integral terms: The 
difference between two Legendre functions appearing in these terms may be globally bounded (in absolute value) as 



|P m (cos0) - P_i(cos0)| < C [n(l - I)]" 1 / 2 , 



(A17) 



where Co is a number independent of I and 9 [see Eq. 8.838 of Ref. |27|]. For any 9 and 6o, the integrals of Eq. ( A16| ) 
are thus each bounded (in absolute value) by C\{9, 6q) x l^ 1 / 2 , and we find 



E SG U 



< C 2 (9,9 a )r 1/2 [|P m (cos0)| + |P( 



cos( 



(A18) 
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as I — > oo (the above coefficients C% and C2 are /-independent). Thus, the standard infinite sum over <5Gv e \ vanishes, 
and hence also the tilde-sum: 



(A19) 



With Eqs. (A.14) and (A19), Eq. rtASh is verified 



APPENDIX B: DERIVATION OF INTEGRALS 



In this appendix we obtain the tilde- integrals given in Eqs. ( |99| ) and (100), which are needed in the calculation of 
the regularization parameters. We start with Eq. (p9[). Let I k (X) denote the primitive function of X k J„(A), where A 
is a real variable, J n is the Bessel function of the first kind, of order n, and k, n E N: 



J*(A) = / A k J„(A)dA. 



Let also 1^ stand for the definite integral 



r>A poo 



I k ee limA^oo / {X'Y J n {\')d\' ee / X k J n (X)d\ 



(Bl) 



(B2) 



Consider first the case k = 0. The standard integral J °° J n (X)dX — 1 is well defined and finite. Thus, in this case, 
the tilde-integration in Eq. (p2) may be replaced with a standard integration, yielding 



e=° = l, Vn>0. 



(B3) 



Consider next the case k > 0. Writing in Eq. @ A fc J„(A) = A fc_ ™ _1 [A" +1 J„(A)] = A fe_ ™ _1 [A™ +1 J„+i(A)] ' [where 
use is made of Eq. (pq) and a prime denotes <i/dA], and integrating by parts, we arrive at the recursive formula 



I k n (X) = X k J n+1 (X) -(k-n- l)Ci(A). 
If < k < n, then by k successive applications of this recursive formula we obtain 



3=0 



( n _fc_l + 2 j)!! 1 (n + k-l)\l 

(n-k-l)V A J «+ 1 +^ A )J + („-fc-l)!! Z "+ fe(A) ' 



Now, the Bessel functions J n (A) admit the asymptotic form 

J„(A 00) ~ (2/vrA) 1/2 cos (A - n7r/2 - tt/4) 



(B4) 



(B5) 



(B6) 



(see, e.g., Eq. 8.451-1 in |27f]). Therefore, each of the k terms in the sum over j in Eq. (B5) diverges at large A as some 



positive (half-integer) power of A times oscillations with respect to A. Clearly, all such terms are eliminated when the 
tilde-limit A — > 00 is taken. Also, all of these terms vanish at A = 0. We are thus left with 



j*_ (" + fc-l)!! jt, _ (n + fc-l)!! {QrQ<k<n 
ln ~ (n-k- 1)!! In+k ~ (n-k- 1)!!' iOT ° < ~ 



(B7) 



where the last equality is due to Eq. (B3) 



If k > n and the difference k — n is even, then by p = (k — n)/2 applications of the recursive formula ( B4 ) we obtain 



/*(A) = ^(-iy 
3=0 



(fc-n-l)!! xk 



-J—X k -> J n+1+] {X) + (-ink n 1)H/JKS^(A). 



(ifc-n-1 -2j')! 

Again, all terms in the sum over j vanish at the tilde- limit A — > 00 and at A = 0, leading to 



(B8) 



2G 



(-lHfc-n-l)!!jf£g/ 2 2 (A) = 

{-l) {k - n)/2 {k -n- l)!!(fc + n - 1)!!, for even k - n > 0, 



(B9) 



where the last equality is due to Eq. (B7). 

The situation is differe nt in case k > n and the difference k — n is odd. Then, following q = (k — n—l)/2 applications 
of the recursive formula ( B4 ) one obtains 

i k nW = B-D j (k -n-i-2m xk ~ iJ ^ x) + ( ~ 1)g(fc - n - = 



J=0 



(k-n-1- 2j) 
(k-n- 1)!! 



(A- 



1 - 2j) 



-A fe J J„ +1+i (A), 



(BIO) 



with no residual integral [notice that when Eq. (B4) is applied with the upper index of J„ greater by 1 than its lower 
index, then the second term on the RHS of this recursive formula vanishes]. This leads, when applying the tilde-limit, 
to 



0, for odd k — n > 0. 



(Bll) 



The above results, Eqs. (|BJ), flBTj), flBJ), and ( |Bll| ), are summarized by Eq. (|99|) in Sec. [VII. 
We further need now to calculate the integral given in Eq. ( 100 ). Integrating by parts, we express the primitive 
function of the integrand, A fc J n (A) In A, for odd k — n > 0, as 



jfe(iog) _ / A fc J„(A)lnAdA = ^(A)lnA- / [Z*(A)/A] dX. 



(B12) 



By virtue of Eq. ( B1C ), the surface term on the RHS here is dominated at small A by oc \i k + n + 1 )/ 2 J( fc+n+1 w 2 (A) In A oc 
yk+n+i \ n x, and therefore it vanishes at the limit A — > 0. However, this surface term diverges at the tilde-limit A — > oo 
(as A' c_1 / 2 lnA times oscillations with respect to A). This can be avoided by slightly modifying the definition of a 
function's tilde-limit, by allowing the quantities in Eq. (A6) to also admit the form Bjf* = cijl bj InZ cos(ajl + /3j) 
(with ctj ^ for all j). It can be shown that all features and results discussed in Appendix \A\ concerning the tilde- limit 
remain valid also under this wider definition. With the revised definition of the tilde- limit, the surface term in Eq. 
(B12) vanishes at the tilde-limit A — * oo as well as for A — > 0. Denoting i^ 1 ° s ' ) = limA^oo-^ los \ we then have 



jfe(log) = 



(fe-n-l)/2 



\ [i k n {\)i\\d\ = - Yl 



(k-n- r 



(k 



(B13) 



-l-2j)!! n+1 +^ 

n is an odd number, the difference between 



where we have substituted for i^(A) from Eq. (B10). Provided that k 
the upper and lower indices of i„^7_- is also an odd number, k — n — 2 — 2j. Therefore, by virtue of Eq. (Bll), we 
- i ■ = for any j satisfying k — j — 1 > n + 1 + j, i.e., j < (k — n — 2)/2. We find that of all the terms 



have = 

summed up in Eq. ( B13[ ) , the only nonvanishing one is the one with j = (k ■ 



l)/2. Hence, 



= -(-l) 1 



(fc-„_l)/ 2 , fe _ n _ 1 V| J[k+n-l)/2 = ( _ u (k-n+l)/2 (k 



(fe+n+l)/2 



= (-l) 1 



l)!!(fc + n- 1)!!, 



(B14) 



where, in the last equality, the value of I 



(fc+n+l)/2 



has been inferred from Eq. (B7). This proves Eq. (100). 
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